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Abstract 

The Isolated Horizon formalism, together with a simple phenomenological 
model for colored black holes was recently used to predict a formula for the 
ADM mass of the solitons of the EYM system in terms of horizon properties 
of black holes for all values of the horizon area. In this note, this formula is 
tested numerically -up to a large value of the area- for spherically symmetric 
solutions and shown to yield the known masses of the solitons. 
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Historically, the use of simple models has played an important role to gain intuition 
about the behavior of complex gravitational systems. For example, the intuitive notion 
of "what can be radiated will be radiated" after a gravitational collapse gave way to the 
formulation of the "no-hair conjecture" for stationary solutions Similarly, the fact that 
there exist no solitonic (i.e. stationary, regular) solutions to the Yang-Mills equations on 
Minkowski space-time was in part responsible for the general belief that solitons were also 
absent in the presence of gravity. However, as these examples illustrate, the intuitive models 
that one constructs are of a limited validity once one considers more general situations. For 
instance, it is now known that the Einstein- Yang-Mills system violates both of the early 
"conclusions", that is, there exist both solitons and colored black holes in that 
theory. More recently, the use of dimensional analysis has helped to develop intuition about 
the existence of static (regular and BH type) solutions in theories with different matter 
content 0. The basic idea is that whenever a dimension-full constant can be constructed 
from the fundamental coupling constants of the theory, then there exist non-trivial static 
solutions. For example, in the EYM system the gravitational constant G and the gauge 
coupling constant g give raise to a mass parameter (or a length scale) that governs the scale 
of the mentioned solutions. One of the facts that can be heuristically explained from this 
viewpoint is the emergence of new phenomena once more fields (and thus more coupling 
constants) are introduced. For instance, in the passage from EYM to EYM-Higgs, a new 
mass parameter can be constructed from the expectation value r] of the Higgs field. This 
indicates, for instance, that the hair of the hairy black holes is always short and that those 
BH have an upper bound on their area [0. 

During the past year, the introduction of the Isolated Horizon (IH) formalism has 



allowed to gain insight into the static sector of theories admitting "hair" |TT|. First, it was 
found that the Horizon Mass of the black hole, constructed purely of quasi-local quantities, 
is related in a simple way to the ADM masses of both colored black holes and solitons of 
the theory [^. Second, a simple model of colored black holes as bound states of regular 
black holes and solitons emerged and provided heuristic explanations for the behavior of 
horizon quantities of colored black holes Finally, the combination of the Mass formula, 



together with the fact that in EYMH theory, different "branches" of static solutions merge, 
has lead to a formula for the difference of soliton masses in terms of black hole quantities 
I3| , p!2[] . In this article, we explore further the consequences of the IH formalism in the 



static sector of the theory. In particular, we explore a formula for the ADM mass of the n 
Bartnik-McKinnon soliton in terms of horizon quantities of the related black holes, 

M^S = ^ f {I - Mr))dr (1) 

where (r) = 2r«;(„) (r) and (r) is the surface gravity of the n colored black hole of 
area a a = 47rr^. This formula is striking because it allows to compute masses of solitons 
from surface gravity of black holes. It should be noted that even when this formula refers to 
the EYM system, it was predicted from considerations in the EYMH system ||1 2|| . 

Two basic ingredients are needed for the validity of this formula: i) That the static sector 
of the EYMH system, in the limit when the vacuum expectation value rj approaches zero 
(t) i-^ 0) reduces to the static sector of the EYM system. Thus, the solutions of the EYM 
system would correspond to limiting solutions of EYMH equations; ii) The surface gravity 
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of colored black holes should have certain limiting behavior as the horizon area grows (in 
order to have convergence of the integral in (|l|)). In this note we shall prove i) and, by 
numerical integration, we show evidence for the mass formula (|1|), and thus indirectly of the 
limiting behavior of surface gravity. Our numerical integration is restricted to the spherically 
symmetric sector and for the n = 1,2 colored black holes, but one expects the main result 
to be valid in general, non-spherical solutions and for all values of n. 

Let us begin by recalling the Isolated Horizons (IH) formalism, whose applications range 
from the extraction of physical quantities in numerical relativity to quantum entropy calcu- 
lations 1^. The basic idea is to consider space-times with an interior boundary (to represent 
the horizon), satisfying quasi-local boundary conditions ensuring that the horizon remains 
'isolated'. Although the boundary conditions are motivated by geometric considerations, 
they lead to a well defined action principle and Hamiltonian framework. Furthermore, the 
boundary conditions imply that certain 'quasi-local charges', defined at the horizon, remain 
constant 'in time', and can thus be regarded as the analogous of the global charges defined 
at infinity in the asymptotically flat context. The isolated horizon Hamiltonian framework 
enables one to define the Horizon Mass Ma, as function of the 'horizon charges'. 

In the Einstein-Maxwell and Einstein-Maxwell-Dilaton systems considered originally , 
the horizon mass satisfies a Smarr-type formula and a generalized first law in terms of 
quantities defined exclusively at the horizon (i.e. without any reference to infinity). The 
introduction of non-linear matter fields like the Yang-Mills field brings unexpected subtleties 



to the formalism [jrT|. However, one can still define a Horizon Mass, and furthermore, this 
Horizon Mass satisfies a first law. 

An isolated horizon is a non-expanding null surface generated by a (null) vector field 1°'. 
The IH boundary conditions imply that the acceleration k of {l°''VJ'' = kI^) is constant 
on the horizon A. However, the precise value it takes on each point of phase space is 
not determined a priori. On the other hand, it is known that for each vector field on 
space-time, the induced vector field Xt^ on phase space is Hamiltonian if and only if there 
exists a function Et^ such that 6Et^ = Q{S,Xt^), for any vector field 6 on phase space. This 
condition can be re-written as 0], 6Et^ = Soa + work terms. Thus, the first law arises 
as a necessary and sufficient condition for the consistency of the Hamiltonian formulation; 
the allowed vector fields t"" will be those for which the first law holds. Note that there are as 
many 'first laws' as allowed vector fields = t"" on the horizon. However, one would like to 
have a physical first law, where the Hamiltonian Ef^ be identified with the 'physical mass' Ma 
of the horizon. This amounts to finding the 'right k\ This 'normalization problem' can be 
easily overcome in the EM system In this case, one chooses the function k = /c(aA, Qa) 
as the corresponding function for the static solution with charges (aajQa)- However, for 
the EYM system, this procedure is not as straightforward. A consistent viewpoint is to 
abandon the notion of a globally defined horizon mass on Phase Space, and to define^ for 
each value of canonical normalization t^ that yields the Horizon Mass for 



the Ho branch [^. The horizon mass takes the form. 



Mi"°)(i?A) = ^/^ VK)(r)dr, (2) 
along the Hq branch with Ra the horizon radius. 
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Finally, for any n one can relate the horizon mass M^"^ to the ADM mass of static black 
holes. Recall first that general Hamiltonian considerations imply that the total Hamiltonian, 
consisting of a term at infinity, the ADM mass, and a term at the horizon, the Horizon Mass, 
is constant on every connected component of static solutions (provided the evolution vector 
field t"" agrees with the static Killing field everywhere on this connected component) P,p!0 



In the Einstein- Yang-Mills case, since the Hamiltonian is constant on any n-branch, we can 
evaluate it at the solution with zero horizon area. This is just the soliton, for which the 
horizon area oa, and the horizon mass Ma, vanish. Hence we have that if^*"'"-) = Mg"j\ 
Thus, we conclude: 

M^} = M^i]m - Mi") (3) 



on the entire nth branch |]TT],|TO|. Thus, the ADM mass contains two contributions, one 
attributed to the black hole horizon and the other to the outside 'hair', captured by the 
'solitonic residue'. The formula (^, together with some energetic considerations, lead to 
the model of a colored black hole as a bound state of an ordinary, 'bare', black hole and 
a 'solitonic residue', where the ADM mass of the colored black hole of radius -Ra is given 
by the ADM mass of the soliton plus the horizon mass of the 'bare' black hole plus the 



binding energy |T2|. Simple considerations about the behavior of the ADM masses of the 
colored black holes and the solitons, together with some expectations of this model (such 
as a negative binding energy) give raise to prediction of the behavior of horizon parameters 



1^ . Among the predictions, we have that /3(n)(-RA); as a function of -Ra and n, is a positive 
function, bounded above by /?(o)(-Ra) = 1- Besides, the curves /3(„)(r), as functions of r 
intersect the r = axis at distinct points between and 1, never intersect, and have the 
property that the higher n is, the lower the curve. Finally, the curves, for large value of 
their argument, become asymptotically tangential to the curve /?(o)(-Ra) = 1- One of the 
features of these solutions is that there is no limit for the size of the black hole. That is, if 
we plot the ADM mass of the BH as function of the radius /2a we get an infinite number 
of curves, each of them intersecting the R/^ = line at the value of the soliton mass, and 
never intersecting each other. 

When one introduces extra fields, as is the case of a Higgs field in EYM-Higgs, then 
there is an extra dimension-full parameter r] that brings a new mass scale into the problem. 
This, in turn, has as a consequence that there is a maximum value for the horizon radius 
of the black holes |l7|JT^. This also indicates that in the Madm vs. -Ra plot, the curves 
corresponding to two families of black holes intersect at a point. Now, recall that the 
horizon Mass was obtained by integrating M^°\R/^) = ^ J^'^ (^{no){^) d?", along the curve 
corresponding to the nth family. If now two such curves, say the n-th and the n + 1-th, 
intersect, one can use this formula together with Eq. (|^), to conclude that, 

Mir^^-MfJ = ^//3(r)dr (4) 
where the closed counter integral is performed by first moving along the n-th branch up to 



the crossing point and then returning along the n + 1-th branch back to -Ra = ||T3|,|T2 
This concludes the summary of the IH formalism needed for this note. 

Now let us consider the static sector of a family of EYMH theories, parametrized by rj. 
As T] decreases, the second mass parameter of the theory becomes smaller and the upper 
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bound on the area of the hairy black holes increases |^. Thus we expect that the point 
of intersection of the n and n + 1 branches would move towards larger values of i? as 77 
decreases. Moreover in the limit i] ^ we would expect this intersection point to move 
towards infinity leading to a situation where the different branches do not intersect, as is in 
fact observed to happen in pure EYM theory. We might also argue that as r] determines 
the natural scale for the Higgs field in static situations, we expect that as 77 — * the static 
sector of EYMH theory would have vanishing Higgs field and thus correspond to the static 
sector of pure EYM theory. We will support this argument by explicitly proving that for 
the case 77 = the static, purely magnetic solutions of EYMH theory have vanishing Higgs 
field. The proof is a simple generalization of a Bekenstein no hair theorem [ p!5[] . 

Consider an EYMH theory with scalar field $ with potential V{^) = A($*<l> - r/^)^, 
where A , and rj are constants, and is the Hermitian conjugate of $. Consider a static 
black hole solution with time-like Killing field The equation of motion for the scalar field 
is: 

dV 

where Da stands for the gauge and metric covariant derivative Da = Va — ie.Al^T\ with 
Va is the metric compatible derivative operator, stand for the gauge fields, for the 
Lie algebra generators, and e is the gauge coupling constant. For non-extremal black holes 
without loss of generality we can consider that the space-time has a bifurcate Killing Horizon 
with bifurcation surface 5" [|14| . Let t be the Killing parameter and consider M the region of 



space-time bounded by Si and S2 hyper-surfaces of constant t, S and asymptotic infinity. 
We multiply Eq. by $* and integrate over M: 

r 8V 
0= dH'^*D^Da<i>-<i>*—-)V^ (6) 



M 9$* ■ 

= / di:''^*Da^- [ d'^x^^(Da^*Dh^g'''' 

JdM JM 

Consider the case where 77 = 0. The boundary integral consists of four terms: The integrals 
over El and S2 that are equal in magnitude and opposite in sign due to the time translation 
invariance and the vanishing of the integral at infinity because of the fall-off conditions on 
$ required from asymptotic flatness. Finally, the term associated with S which does not 
contribute since the integral is over a lower dimensional manifold. 

Next we write the inverse metric g"''' = -A^^^^^ + where iV"^ is the norm of the 
Killing field and /i"* is the inverse of spatial metric on S. Using the fact that (^"Va$ 
vanishes as long as the solution is static, and that Aq = ^""Al^ = corresponding to the 
purely magnetic sector, we obtain: 

r 3V 

/ d^x^g{Da^*D,^h''' + ^*—-) = Q (7) 

The first term in the integrand is positive semi-definite in general and the second term 
becomes positive semi-definite when 77 = 0. Thus in this situation the only possibility is 
Da^ = and = which for 77 = implies $ = 0. This is what we wanted to show. 
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TABLE I. R is the upper limit of the integration interval. I(n = 1), I(n = 2), 
I{n = 2) — I{n = 1) are the values of the integral for one node, two nodes, and the difference 
of both, respectively. The complete output of the simulation is reported even when the numerical 
errors make some of the last figures irrelevant. 
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Thus, by considering the limit of the static and purely magnetic sector of EYMH theory in 
the limit rj ^ Owe expect to obtain the static and purely magnetic sector of the EYM theory. 
Furthermore, by considering Eq. in this limit we obtain an equivalent equation (0), now 
for the EYM case. Let us now consider Minkowski space-time as the zeroth (degenerate) 
soliton, where the Schwarzschild branch (the 0**^ branch) begins, we obtain a formula for the 
mass of the nth soliton in EYM theory: 

Mi:^ = ^fp{r)dr, (8) 

where the integral is performed first along the Schwarzschild branch to r = cxo and returning 
to r = along the n}^ branch. This is precisely Eq. (|I]). 

Now we want to establish the validity of formula (0) for EYM solitons. For that, we show 
a comparison between the numerical evaluation of the masses of the EYM solitons and the 
corresponding evaluation of the left hand side of Eq. (|ID for the n = 1,2 colored black holes. 
To do this, we calculate numerically the functions I3{n){r) for n = 1, 2 as a function of the 
radius of the horizon i?A by using a fourth-order Runge-Kutta algorithm for the computation 
of the configurations of black holes with n = 1,2 nodes maintaining the truncation error 
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below of the accuracy 10~^ for the metric and the Yang-Mills field of each configuration. At 
the same time, a fifth-order Runge-Kutta algorithm was used for estimating the error in the 
values of the function /3(„)(r). Then we calculate the value of the integral (|I]) for different 
values (sufficiently large) of the upper limit of integration R and the respective accumulated 
error associated with it coming of the one associated with the values of the function /3(„)(r). 
The accumulated errors associated with the integral (|l|), calculated as the difference between 
the fourth and fifth order algorithms, for n = 1,2 and superior limit of integration R = 150, 
are 5 x 10~^ and 9 x 10"'* respectively (the accumulated errors for R > 150 differ of these 
values by a quantity of the order 10~^° and then have practically converged to these). The 
black hole radius is in the standard units given by the EYM coupling constant (where R = 1 
corresponds to the natural length scale given by the theory). Table H shows the values of 
the integral (0) for n = 1,2 and the difference between both as a function of the upper limit 
of integration R, using the fifth-order algorithm. It is important to notice from the Table 
that the difference of the integrals, which by equations (|I|) and (^ should approximate the 
difference in soliton masses, has indeed a faster convergence than the integrals by themselves. 
We have restricted the upper limit in the integral to R = 150 given that the difference in 
masses had practically converged to the reported value. 

From the values given in Table |, it is clear that. For R sufficiently large,, the integrals for 
n = 1,2 (and the difference of both) are approaching to the corresponding soliton masses, 
0.82864698216 and 0.97134549426 for n = 1,2 respectively (and the difference of them. 



0.142698512), as reported in p 



To conclude, within the uncertainty of our numerical calculation, the soliton mass formula 
(|I|) provides the ADM mass of the EYM system. 
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